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Abstract 



Oh 



■ Using an effective strongly coupled lattice QCD Hamiltonian 
CN . and Wilson fermions we calculate the equation of state for cold 

and dense quark matter by constructing an ansatz which exactly 
• diagonalizes the Hamiltonian to second order in field operators for 

all densities. This ansatz obeys the free lattice Dirac equation with 
a chemical potential term and a mass term which is interpreted as 
the dynamical quark mass. We find that the order of chiral phase 
transition depends on the values of input parameters. In the phase 
with spontaneously broken chiral symmetry the quark Fermi sea 
I has negative pressure indicating its mechanical instability. This 

■ result is in qualitative agreement with those obtained using contin- 
uum field theory models with four-point interactions. 

• 1—1 - 

X 

1 Introduction 

Simulating finite density QCD is one of the outstanding problems in lat- 
tice gauge theory Because of the sign problem no reliable numerical 
simulations of finite density QCD with three colors exist even in the strong 
coupling limit || .[] This is a rather frustrating situation in view of the cur- 
rent intense interest in finite density QCD fueled by the phenomenology of 
heavy ion collisions, neutron stars, early universe and color superconduc- 
tivity. Therefore even a qualitative description of strongly coupled QCD 
at finite density using field theoretical methods is quite welcome. 

1 Simulation of two-color four-flavour QCD at finite quark number density has re- 
cently been reported in H . 
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Strongly coupled QCD at finite quark chemical potential fi has pre- 
viously been studied analytically both in the Euclidean |4|, |], U and in 
the Hamiltonian 0, [|] formulations. Except for J7J, the consensus is that 
at zero or low temperatures strong coupling QCD undergoes a first order 
chiral phase transition from the broken symmetry phase below a critical 
chemical potential fic to a chirally symmetric phase above /ic- In this 
work we calculate the equation of state of strongly coupled QCD at finite 
\i using Wilson fermions and show that in the phase with spontaneously 
broken chiral symmetry the quark Fermi sea has negative pressure indi- 
cating its mechanical instability. This result is in qualitative agreement 
with those obtained using continuum field theory models with four-point 
interactions R PHI . 



2 Effective Hamiltonian and the ansatz for finite \i 

We begin by defining our effective Hamiltonian for strongly coupled QCD 



which was first derived by Smit |TTJ. Henceforth we adopt the notation 



of Smit set the lattice spacing to unity (a = 1) and work in momen- 
tum space. The charge conjugation symmetric form of Smit's effective 
Hamiltonian with a chemical potential /zq in momentum space is 



H, 
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where (Ej) = — i (707; — ^7o) with the Wilson parameter r taking on val- 
ues between and 1. In the above Hamiltonian color, flavor and Dirac 
indices are denoted by (a, b), (a, (3) and (a, u, 7, 5), respectively, and sum- 
mation convention is implied. Three of the four parameters in this theory 
are the Wilson parameter r, the current quark mass M and the effec- 
tive coupling constant K which behaves as 1/g 2 with g being the QCD 
coupling constant. The fourth parameter /x is in general not the total 
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chemical potential /x to t of the interacting many body system. As shown 
below the interaction will induce a contribution to /x tot which is momentum 
dependent. We shall therefore refer /x as the "bare" chemical potential. 

Our strategy is to construct an ansatz which obeys the equation of 
motion corresponding to if e ff at finite /x and use it to calculate the equation 
of state of cold and dense quark matter. In free space this ansatz has the 
same structure as the free lattice Dirac field and obeys the free lattice 
Dirac equation with a mass term which is interpreted as the dynamical 
quark mass |12|] . Temporarily dropping color and flavor indices, this ansatz 
is given by 

* v (t,p) = b(p)Up)e- Mf)t + dt(_p)^(_p) e +^ (2) 

The annihilation operators for particles, b, and anti-particles, d, annihilate 
an interacting vacuum state and obey the free fermion anti-commutation 
relations. The properties of the spinors £ and i] are given in |12| . The 



equation of motion for a free lattice Dirac field fixes the excitation energy 
oj{p) to be 

w (p) = (E sin2 (P • ni) + M 2 (p)j (3) 

where M(p) plays the role of the dynamical quark mass or, equivalently, 
the chiral gap. 

We now introduce temperature T and \i simulatneously into our for- 
malism and take the T — > limit to construct an ansatz for finite /x. This 
is accomplished by subjecting the b and d operators in Eq. @ to a gen- 
eralized thermal Bogoliubov transformation as in thermal field dynamics 

H 

b{p) = a p B(p)-(3 p B\-p) (4) 
d(p) = lp D{p) - 5 p D\-p) (5) 

The thermal field operators B and B^ annihilate a quasi-particle and 
create a quasi-hole at finite T and /x, respectively, while D and D* are the 
annihilation operator for a quasi-anti-particle and creation opertor for a 
quasi-anti-hole, respectively. 

The thermal annihilation operators annihilate the interacting thermal 
vacuum state | Q{T,p)) for each T and /x. 

B@)\g(T,n))=B(?)\g(T,n)) = (6) 
D(f)\g(T,n))=D(f)\g(T,n)) = (7) 
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The thermal doubling of the Hilbert space accompanying the thermal 
Bogoliubov transformation is implicit in Eqs. (§) and (0) where a vacuum 
state which is annihilated by thermal operators B, B, D and D is defined. 
Since we shall be working only in the space of quantum field operators it 
is not necessary to specify the structure of \Q(T,fi)). 

The thermal operators satisfy the Fermion anti-commutation relations 
just as the b and d operators in the free space ansatz. The coefficients of 
the transformation are a p = Jl — n~ , (3 P = \Jn^-, 7 P = \Jl — n+ and 

5 P = ynf, where = [e ( ^ ±/ ^/( fcsT ) + are the Fermi distribution 
functions for particles and anti-particles. They are chosen so that the 
total particle number densities are given by 



n, 



n. 



{G{T,ii)\tf(p)b(p)\G{T,ii)) 
{Q{T^)\d\p)d{p)\g{T^)) 



(9) 



Hence in this approach temperature and chemical potential are introduced 
simultaneously through the coefficients of the thermal Bogoliubov trans- 
formation and are treated on an equal footing. We stress that the chemical 
potential appearing in the Fermi distribution functions is the total chem- 
cial potential of the interacting many body system which is not necessarily 
equal to the bare chemical potential /xo appearing in the effective Hamil- 
tonian in Eq. ([]]). 

We demand that our ansatz satisfies the equation of motion corre- 
sponding to the free lattice Dirac Hamiltonian with a chemical potential 
term. This Hamiltonian is given by 
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As in [0 the mass M(p ) will be identified with the chiral gap while /itot 
is the total chemical potential mentioned above. In the T — > limit the 
excitations of anti-holes are supressed which amount to setting j p = 1 
and 5 P = in Eq. (El). Thus our ansatz at finite /i to t is given by 



a p B{p) - l3 p B\-p\ 



+ l p D\-p)r] u (-p)e 



+i[u)(p)+fi tot ]t 



(11) 
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The spinors £ and r\ obey the same properties as in free space and the 
excitation energy u(p) has the same form as in Eq. @. 

Two quantities of interest that can be calculated in a straightforward 
manner using Eq. (|TT| ) are the quark number density ( 1 $ r 7o\E r ) and the 
chiral condensate (tyty. They are found to be 

(* 7 o*> = £/3 p 2 (12) 
p 

and 

We see immediately from Eq. flllf ) that the chiral condensate is propor- 
tional to the dynamical quark mass and therefore can clearly be identified 
as being the order parameter for the chiral phase transition at finite den- 
sity. 



3 Application of the Equation of Motion 

3.1 The Gap Equation and the Induced Chemical Potential 

We proceed by using the ansatz shown in Eq. ( |TTD in our effective Hamil- 
tonian to determine the equation of motion. This is accomplished by 
exploiting the fact that our ansatz satisfies the equation of motion cor- 
responding to the free Hamiltonian at finite fi given in Eq. flUf). We 
therefore have the relation 



(^ aa ) u (t,q),H' 



(^aa) u (t,q),H eS 



(14) 



where the symbol : : denotes normal ordering with respect to the thermal 



vacuum \ Q(T = 0, //)). Evaluating both sides of Eq. fli~4]) and equating 
terms which are linear in the field operators we obtain 



sin(g • fii)(-foli)u8 + M(q ) (70)^5 - /itot<W (^aa) s (t, q) = 

1 

MfiiWiU + B(q)(j )»8 + C(q)dJ q. 



(15) 
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The momentum dependent coefficients B and C in Eq. ( jT5|) determine the 
gap equation and the total chemical potential, respectively. 

M(q) = B(q) 

+ ^E (* - /?) ( 8r2 cos ^' ^) cos ^" ^ 



- i(l + r 2 ) cos(p + <f ) ■ hi ) ( 1 G ) 



c p,i 

The structure of the gap equation Eq. (|16D is very similar to that 
in free space (/5 2 = 0) found in [I2|. The dynamical quark mass is a 



constant to lowest order in N c but becomes momentum dependent once 
1/N C correction is taken into account. We see from Eq. (|i~7l) that the 
total chemical potential is a sum of the bare chemical potential /xq and a 
momentum dependent contribution generated by the interaction. Because 
of this induced chemical potential Eqs. (|I6| ) and ( |TTD are coupled and must 
be solved self-consistently. It should be noted that this shifting of the bare 
chemical potential by the interaction is not a new effect. In the Nambu- 
Jona-Lasinio model [|14) within the Hartree-Fock approximation and at 
finite T and /i, the interaction induces a contribution to the total chemcial 
potential which is proportional to the number density ||15|| . 

Solutions for the coupled equations Eqs. ([TT]) and ( |T7D to lowest order 
in Nq are shown in Figure [1] where we present the dynamical quark masses 
as functions of the total chemical potential for various combinations of the 
Wilson parameter r and the coupling constant K with a vanishing cur- 
rent quark mass. As shown in Figure la we find first order chiral phase 
transitions when r = (i.e. with naive fermions) with critical chemical po- 
tentials of [j,c ~ 1-04 and 1.33 for K = 0.8 and 0.9, respectively. As far as 
the order of the phase transition is concerned, our results agree with those 
obtained using the Euclidean formulation with Kogut-Susskind fermions 
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Figure 1: a) Dynamical quark mass versus the total chemical potential to lowest 
order in N c for K = 0.8 (solid line) and 1.0 (dashed line) with r = 0. A 
first order chiral phase transition takes places in both cases with a critical 
chemical potential of nc = 1-042 and 1.329 for K = 0.8 and 1.0, respectively, 
b) Dynamical quark mass versus the total chemical potential to lowest order 
in N c for K = 0.8 (solid line) and 0.9 (dashed line) with r = 0.25. A second 
order chiral phase transition takes place for K = 0.8 with a critical chemical 
potential of fic = 0.836 while the transition is first order for K = 0.9 with \ic 
= 1.086. The current quark mass is set to zero (Mq = 0) in all cases. 
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@> H? an d also with a recent result obtained with naive fermions (r = 
0) in a Hamiltonian formulation ||. However we disagree with the result 
of Le Yaouanc et al. |7|] who find a second order chiral phase transition 
using the same effective Hamiltonian as in this work and naive fermions. 

For r = 0.25 we see from Figure lb that the phase transition can be 
either first or second order depending on the value of the effective coupling 
constant K. When K = 0.8 we find a second order phase transition with 
He ~ 0.84 while if the coupling constant is increased to 0.9 the phase 
transition becomes first order with a larger critical chemical potential of 
fic ~ 1-09. It would be interesting to investigate whether the approaches 
adapted in f7| and [|J to study strongly coupled QCD in the Hamiltonian 
lattice formulation using naive fermions would result in a similar behaviour 
of the order of the phase transition when extended to r > 0. Finally, we 
note that for both r = and 0.25 the value of the critical chemical potential 
increases as K is increased. The same behaviour has been observed in 0. 



3.2 Off— diagonal Hamiltonian and the Vacuum Energy Density 

The off-diagonal Hamiltonian which is bilinear in field creation and anni- 
hilation operators may be written as 

tfofrl 0(0, //)) = £ F(q)Bl a (q)Dl a (-q) | 0(0, //)) (18) 



where F{q) is a scalar function of q. We see from Eq. ([IB] ) that the 
elementary excitation of the effective Hamiltonian involves color singlet 
(quasi-) quark-anti-quark excitations coupled to zero total three momen- 
tum. With the use of the equation of motion and properties of r] spinor 
[ PI it can be shown analytically that H oS \ Q(0,fi)) = for all fi. There- 
fore our ansatz exactly diagonalizes the effective Hamiltonian to second 
order in field operators for all densities. 

Having verified that our ansatz diagonalizes the second order Hamil- 
tonian we are in a position to evaluate the vacuum energy density. Using 
Eq. fllj) we find 

i(0(o,/i)|if eff |0(O,/i)) = 



al 



!^l+r 2 )+^ + ^M 
2 u(p) 
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+ (l + /5> (19) 



From Eq. (f03) we can evaluate the difference of the vacuum energy densi- 
ties in the Wigner-Weyl {M(q) = 0) and Nambu-Goldstone (M(q) ^ 0) 
phases of the theory. We find 

ae = ^(g\H cS \ g)\ Am=0 -^(G\Hes\ G)\ M m?o > (20) 

and therefore the true ground state of our interacting many body system 
is in the phase with broken chiral symmetry. 



4 Equation of State 

Having solved the gap equation for the dynamical quark mass we proceed 
to determine the equation of state of cold and dense quark matter. We 
present results obtained with a vanishing current quark mass and Wilson 
fermions for r = 0.25 to leading order in N c since in this case we observe 
both first and second order chiral phase transitions. The pressure of the 
cold and dense quark matter is obtained by numerically evaluating the 
thermodynamic potential density. In Figure [| we plot the pressure as 
functions of quark number density for K = 0.8 and 0.9. In both cases 
we find the pressure to be negative in the broken symmetry phase and 
therefore the quark Fermi sea is mechanically unstable with an imaginary 
speed of sound. 

Our conclusion regarding the quark matter stability at finite density is 
consistent with results obtained using the Nambu-Jona-Lasinio model [[J 
and the effective instanton induced 't Hooft interaction model [[Hj . In both 
M and |TD| mean field calculations show that cold and dense quark matter 



is unstable in the phase with spontaneously broken chiral symmetry and 
lead the authors to speculate the formation of quark droplets reminiscent 
of the MIT bag model. 

5 Conclusion 



In this work we studied strongly coupled QCD in the Hamiltonian lat- 
tice formalism at finite density using Wilson fermions. Starting from an 
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Quark Number Density 



Figure 2: Pressure versus the quark number density to lowest order in N c for 
K = 0.8 (solid line) and 0.9 (dashed line) with M = and r = 0.25. Quark 
number densities at the chiral phase transtion point are n q = 0.013 and 0.025 
for K = 0.8 and 0.9, respectively. 



effective Hamiltonian we constructed and an ansatz which exactly diago- 
nalizes the Hamiltonian to second order in field operators for all densities. 
This ansatz obeys the free lattice Dirac equation with a chemical potential 
term and a mass term which plays the role of the dynamical quark mass. 
This mass and the total chemical potential of the interacting many body 
system were determined to lowest order in N c by solving a coupled set 
of equations obtained from the equation of motion. We find a first order 
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chiral phase transition with naive fermions while the order of the phase 
transition can be either first or second for Wilson fermions. As in earlier 
studies using continuum four-fermion interaction models we find that the 
quark Fermi sea is mechanically unstable in the phase with broken chiral 
symmetry. Our results for the equation of state of cold and dense quark 
matter should certainly be verified in any future numerical simulation of 
finite density QCD. 
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